Random Organization and Plastic Depinning 
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We provide evidence that the general phenomenon of plastic depinning can be described as an 
absorbing phase transition, and shows the same features as the random organization which was 
recently studied in periodically driven particle systems [L. Corte et al., Nature Phys. 4, 420 (2008)]. 
In the plastic flow system, the pinned regime corresponds to the absorbing state and the moving 
state corresponds to the fluctuating state. When an external force is suddenly applied, the system 
eventually organizes into one of these two states with a time scale that diverges as a power law at 
a nonequilibrium transition. We propose a simple experiment to test for this transition in systems 
with random disorder. 

PACS numbers; 64.60.Ht,83.50.-v,62.20.Fk 



In plastic depinning in two dimensional (2D) systems 
with random disorder, which appears in a wide range 
of different systems, particle motion occurs in the form 
of intricate fluctuating channels in which some particles 
are mobile while others remain pinned P, This de- 
pinning process has been studied extensively over the 
years in many systems, including vortices in type-II su- 
perconductors Hi S S Hi vortex motion in Josephson- 
junction arrays [5|, plastic depinning in electron crystals 
on random landscapes @ , charge transport in disordered 
metallic dot arrays 0], and fluid flow on a rough sub- 
strate [1, 0] . Other systems which exhibit similar behav- 
ior include the unjamming and depinning of dislocations 
lol 11 1 and the motion of magnetic domain walls [l^ . 
Complex plastic flow patterns were observed directly at 
the depinning transition in recent simulations and ex- 
periments on a model system of colloids interacting with 
random pinning Evidence that plastic depin- 

ning is a nonequilibrium phase transition includes the 
fact that the velocity force curves scale a.s v (x {F — FcY^ , 
where Fc is the depinning threshold. Many simulations 
of filamentary plastic flow find /?> 1.5 0,13, S E [li, 
but this exponent has not been explained theoretically. 
The true nature of the pinned to plastic flow transition 
is still not fully understood despite the ubiquity of the 
phenomenon. 

The pinned to plastic flow transition can be character- 
ized as a transition from a nonfluctuating (pinned) state 
to a to fluctuating (plastic flow) state. The same type 
of transition from a nonfluctuating to a fluctuating state 
was recently examined by Corte et al. in a system of 
sheared colloids, which exhibits diffusive liquid behavior 
rather than plastic flow in the fluctuating state [13] ■ A 
novel protocol was used to show that there is a nonequi- 
librium continuous phase transition between the fluctuat- 
ing and nonflutuating states. The system was suddenly 
subjected to a shear and the amount of time required 
for the system to organize into either the fluctuating or 
nonfluctuating state diverged as a power law at the tran- 
sition. Since the system is always spatially disordered 



or random, Corte et al. termed the transition "random 
organization." By using the same type of protocol, here 
we show that plastic depinning exhibits exactly the same 
critical behaviour found in the shearing work, strongly 
suggesting that the transitions in the two systems fall 
into the same universality class. 

In Ref. [13], Corte et al. performed a 2D colloidal 
shearing simulation and suddenly applied a periodic 
shearing force to a collection of colloids with short range 
interactions. The colloids experience a random force 
when they collide. The system always starts in a fluctu- 
ating state where the colloids are diffusing, and organizes 
into either a fluctuating or nonfluctuating state after a 
transient time r. This transient time diverges as a power 
law at a well defined shearing threshold. The experi- 
ments of Ref. also exhibit a power law divergence 
of the transient time at the transition; however, the ex- 
ponents are smaller than those found in the simulation 
which may be due to the three dimensional (3D) nature 
of the experiments. From the exponents and the general 
behavior of the shearing system, the transition between 
the fluctuating and nonfluctuating states is most consis- 
tent with an absorbing phase transition in the universal- 
ity class of directed percolation fl^; however, it has been 
proposed that this particular system falls into the class of 
conserved directed percolation (C-DP) since the number 
of colloids is fixed [201 ]. 

Here we use a similar protocol to show that the transi- 
tion from pinned to plastic fiow exhibits the same power 
law diverging transient times as the system organizes 
into either a nonfluctuating pinned state or a fluctuat- 
ing plastic flow state. The exponents of the divergence 
are close to those observed in the colloidal shearing sys- 
tem of Ref. [13] ■ We also examine quantities that were 
not measured in the shearing work and show that there 
is a power law decay in the number of active particles 
at the transition with an exponent that is in agreement 
with that predicted for 2D (C-DP). We note that the 
idea of elastic depinning transitions falling into a class of 
absorbing phase transitions was previously proposed for 
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the depining of elastic lines which fall into the class of ID 
(C-DP) [2l[; however, this is a very different system from 
plastic depinning since there is no tearing. Additionally, 
the exponents for the elastic system arc very different 
from what we find. 

We specifically model a 2D colloidal system with ran- 
dom quenched disorder, where it has been well estab- 
lished in simulations and experiments that a plastic flow 
phase occurs for sufficiently strong disorder [3, fl^ . 
Since the colloid-colloid interactions are relatively short 
ranged, we can perform much larger simulations for much 
longer times near the depinning transition than have pre- 
viously been done for 2D plastic flow systems. Previous 
work for a periodically driven vortex system showed a 
transition from reversible to irreversible flow 22|; how- 
ever, in Ref. 22| the particles were moving in both states 
and there was no pinned state, placing the system in a 
distinctly different regime than that studied here. 

We consider a 2D system of size L x L with pe- 
riodic boundary conditions in the x and y directions. 
The sample contains Nc colloidal particles with density 
ric = Nc/ and the time evolution of colloid i at position 
Ri is governed by the overdamped equation of motion 
■qdHi/dt = Ff + Ff -f Frf 3, where the damping con- 
stant 77 = 1. The colloid-colloid interaction potential has 
a Yukawa form, V{Rij) = {Eo/Rij) exp{—KRij), where 
Rij = |Ri - Rjl, Eq = Z*^/(47reeoao), e is the solvent 
dielectric constant, Z* is the effective charge, 1/k is the 
screening length, and the force F?'^ — — X^^i ^^(-^y )■ 
Lengths are measured in units of ao, assumed to be 
of the order of a micron, forces in units of Fq = 
Eo/ao, and time in units of ryag/i^o- We neglect hy- 
drodynamic interactions since we are in the low vol- 
ume fraction, highly charged, electrophoretically driven 
limit and since the pinned colloids provide the equiv- 
alent of a porous medium [2^. The substrate force 



Fs — — J2k=i ^^pi^ik') arises from Np randomly placed 
pinning sites of density Up = Np/ LP' , radius Vp = 0.2, and 



maximum force Fp, with Vp{R^^, 
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8(rp — i?'^"*), where 8 is the Heaviside step function 



'iFp/2rp){R^i 



and i?. 



'ik 



IR, - R 



(p)i 



is the distance between particle i 
,(p) 



and a pin at position Rj. . The driving force Frf = FdSi 



represents the effect of an applied electric field [14|, and 
here we consider Fd = 0.1. The initial colloid configu- 
ration is prepared by simulated annealing with Fd = Q. 
Starting from a fully ordered state does not change the 
long-time response of the system. We measure the total 
velocity V = Y^f" dRi/dt ■ x. We consider two system 
sizes, L ^ 24 and 48, with Uc = 2.9 and Up — 3.0. For 
L = 24, Nc = 1672 and Np = 1727, while for L = 48, 
Nc — 6688 and Np = 6912, the largest system of this 



type studied to date l^, l3| . In the 2D shearing simula- 
tions of Ref. [is'l, 1000 particles were used to explore the 
diverging time scale, so our system should be sufficiently 
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FIG. 1: Colloid positions (black dots) and trajectories (black 
lines) over a fixed time of 5 x 10"^ simulation time steps in a 
sample with L = 24. (a-d) Fp/Fp = 0.93 after (a) 2.5 x 10^, 
(b) 1.5 X 10*, (c) 1.5 X 10^ and (d) 1.5 x lO'^ simulation time 
steps, showing that the initial motion settles into a steady 
fluctuating state, (e-h) Fp/F^ = 1.05 after (e) 2.5 x 10=*, (f) 
1 X 10'', (g) 5 X 10*, and (h) 1 x 10^ simulation time steps, 
showing that the system settles into a pinned state. 



large to capture the divergence. 

In Ref. [l3|, we showed that a system with quenched 
disorder exhibits a regime of plastic flow as a function 
of Fp, and that there is a well defined depinning thresh- 
old Fc as a function of Fd- In this work, we perform a 
series of simulations for different values of Fp in which 
we suddenly apply a constant Fd- The system exhibits a 
transient motion before settling into a completely pinned 
state or steady moving state. Similar results appear for 
a periodic driving force with a sufficiently long period. 
For fixed Fd, there is a critical pinning strength Fp such 
that for all Fp > Fp the system settles into a pinned 
(absorbing) state, while for Fp < Fp the system stabi- 

0.1, F^ = 0.3470; 



lizes in a fluctuating state. For Fd 
thus, the transition from a pinned to a fluctuating state 
is not simply determined by when the driving force is 
higher then the pinning force, but is instead affected by 
the colloid-colloid interactions. In Fig. [^a-d), we show 
an example of a system with Fp/Fp = 0.93, below the 
threshold. Shortly after the drive is apphed. Fig. [T]Ja) 
indicates that a large portion of the system is in motion. 
After 1.5x 10^ simulation time steps, Fig.[IJb) shows that 
the number of moving colloids decreases, while for longer 
times the number of active colloids settles down to fluc- 
tuate around an average value, as shown in Fig. [H^c) at 
1 X 10^ simulation time steps and in Fig.[ljd) at 1.5 x 10^ 
simulation time steps. For Fp/Fp > 1.0, the number of 
active particles decreases to zero over time and the sys- 
tem reaches the pinned state, as illustrated in Fig. [ij e-h) 
for Fp/Fc = 1.05. 

In Fig. [U^a) we plot V{t), the time trace of the to- 
tal velocity, for a system with L = 48 at F^ = 0.35, 
where the system settles into a fluctuating state, and 
Fp = 0.3, where the system reaches a pinned state. To 
measure the transient time, we use the same procedure as 
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FIG. 2: (a) V{t) for a system with L = 4 
Fp = 0.3 and Fp/F^ = 0.86; lower curve: 



i. Upper curve: 
Fp = 0.35 and 
Fp/ Fp — 1.008. In both cases a transient appears which is 
much longer near Fp/Fp — 1. Smooth lines are guides to the 
eye. (b) The transient time r versus Fp to decay to the steady 
fluctuating state (circles) or to the pinned state (squares) for 
the system in Fig. 2 with L — 48. Solid lines: power law fits 
to the form r oc |Fp - F^^l"", where v = 1.37 ± 0.06. (c) The 
curves from (b) on a log-log scale. 



Rcf. 18] and fit the decaying V{t) curves to the function 
V{t) {V° - V')eiip{-t/T)/t°' + where V° is the 
initial velocity and V is the steady state velocity. This 
functional form approaches a power law at the transition 
when r — > cx), and the value of a only becomes relevant 
very close to this transition. Performing the same fit for 
different quantities, such as the fraction of moving parti- 
cles or the diffusion in the transverse direction, produces 
the same results. In Fig. E^b) we plot the transient time 
r versus Fp for the sample in Fig. ^a.) on both sides of 



the transition. A clear divergence in the transient time 
occurs near Fp ~ 0.3470. On both sides of the transition, 
the T versus Fp curves can be fit with a power law form. 



^ pi I 



(1) 



with v — 1.37 ± 0.06. Power law fits appear as solid lines 
in Fig. IHb), and in Fig. ^c) we plot the same data on 
a log-log scale. Our system sizes are adequate to obtain 
exponents of this type [l^, and our exponent is in good 
agreement with the value v = 1.33 found in the 2D shear 
simulation [T"" 
gave V = 1.1 
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|. The shearing experiments of Ref. 
however, since the experiments were con- 
ducted in 3D, the exponents could be expected to differ 
from the 2D case. For example, in conserved directed per- 
colation (C-DP), the relaxation time exponent w 1.29 
in 2D and « 1.12 in 3D, while in directed percolation 



(DP), 
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1.295 and 



Fig. [3] we plot T versus Fp 



; 1.105, respectively [li 
for a sample with L - 



In 

24, 

1.36 ± 0.05. 
= 0.317. For 

even smaller systems with L < 24, the periodic bound- 



where we find the same behavior with v = 
The value of F^ shifts slightly down to F^^ 
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FIG. 3: r versus Fp for the steady fluctuating state (circles) 
and the pinned state (squares) in a system with L = 24. Solid 
lines: power law fits with v = 1.36 ± 0.06. Inset: the same 
curves on a log-log scale. 
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FIG. 4: The L = 48 system from Fig. [^b). (a) n^, the 
fraction of moving particles, versus time for Fp/Fp — 0.927, 
1.008, 1.05, and 1.2 (from top to bottom). Straight line: a 
power law flt with a = 1/2. Dashed line: a power law fit 
with a = 0.125. (b) dy, the displacements in the y directions, 
versus time for Fp/F^ = 0.927, 1.008, 1.05, and 1.2 (from 
top to bottom). Straight line: a power law fit with a = 1/2. 
Dashed line: a power law fit with a = 0.125. 



ary conditions begin to affect the results since it becomes 
possible to stabilize single channels of moving colloids. 

In absorbing phase transitions, it is expected that very 
close to the transition, the transient states should show 
a power law decay with I{t) oc t^" [3, 24 1. We con- 
sider two quantities: the fraction of moving colloids, 
urait) = e(|R,(t) - R,(i - S)\ - A), with 



S = 1000 and A = 0.0005, and the transverse displace- 
ments, dy{t) = |(R,(t) - R,(< - ^)) • y|. In Fig.Ha) 
we plot the time decay of n„(t) at Fp/F^ = 0.927, 1.008, 
1.05 and 1.2, while in Fig.[4fb) we show the decay of dy{t) 
for the same pinning strengths. In both cases, curves 
close to the transition can be fit to a power law decay 
with a — 1/2, as indicated by the solid lines in Figs. 
|31[a,b). A similar power law decay with a = 1/2 has 
recently been observed for systems exhibiting absorbing 
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phase transitions in 2D 2J| . Due to the anisotropy intro- 
duced by the driving force, it is possible that our system 
could exhibit ID C-DP; this would give a ~ 0.125, and 
we plot a line with this slope in Fig. [H^b) for compar- 
ison. It would be interesting to examine similar power 
law decays in the colloidal shearing simulations and ex- 
periments. We have also measured the order parameter 
exponent /3 = 0.6 ± 0.06. In 2D, DP gives a w 0.451 and 
/3 « 0.583 while C-DP gives a « 0.50 and /3 « 0.64 [S]; 
our exponents are not accurate enough to distinguish be- 
tween these two models. We are unable to perform a 
finite size scaling analysis since systems with Nc > 7000 
take a prohibitively long time to simulate while systems 
with Nc < 1500 suffer from persistent channels of flow 
induced by the periodic boundaries. Finite temperature 
can smooth the transition, but it should remain observ- 
able at sufficiently low temperatures {25t] . 

Our results could be directly tested in experiments 
with colloidal particles [l3|. It is possible to trap col- 
loids optically |i26|] , so it should also be possible to create 
a disordered optical trap array in which the strength of 
the traps can be controlled by the adjusting the laser 
strength. In this way, a constant drive could be applied 
to the system while the trap strength is varied. Addi- 
tionally, it should be straightforward to conduct similar 
experiments for vortices in samples that are in the plastic 
flow regime. For example, experiments performed near 
the peak effect where plastic flow is expected to occur 
have revealed a decayin g vo ltage signal when a driving 
force is suddenly applied [27| . The voltage is proportional 
to the average vortex velocity, indicating that the vortices 
are undergoing transient motion. With this type of sys- 
tem, the effectiveness of the pinning changes sharply with 
the applied magnetic field near the peak effect; thus, it 
should be possible to apply a constant current and de- 
crease or increase the magnetic field in order to observe 
a transition between a regime that decays into a pinned 
state to one that decays to a fluctuating state, and then 
measure the transient times near the transition. 

In summary, we have shown that a 2D plastic depin- 
ning system exhibits behavior that is very similar to that 
of the recently studied random organization phenomenon 
observed in 2D sheared particle systems. Under a sudden 
applied drive, there is a well-defined critical transition 
where the system organizes into either a non-fluctuating 
(absorbing) state or into a fluctuating state. The tran- 
sient time diverges as a power law at the transition with 
exponents that are the same as those found in a system 
subjected to shear. Our results provide evidence that 
plastic depinning falls into the class of absorbing phase 
transitions which include directed or conserved directed 
percolation. We also propose simple experimental tests 
for this transition in colloidal and vortex systems. 

This work was carried out under the auspices of the 
NNSA of the U.S. DoE at LANL under Contract No. 
DE-AC52-06NA25396. 
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